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Einstein Spacetimes with Weak Regularity 

Philippe G. LeFloch 



Abstract. We review recent work on the Einstein equations of general 
relativity when the curvature is defined in a weak sense. Weakly regular 
spacetimes are constructed, in which impulsive gravitational waves, as 
well as shock waves, propagate. 



1. Introduction 

The Einstein field equations relating the Ricci tensor Ric g , the scalar 
curvature R g and the energy-momentum tensor T of a spacetime (M, g) 

Ry 

2 

can be imposed in a weak sense, provided the metric g belongs to the Sobolev 
space H 1 of square-integrable functions (and satisfies a suitable L°° condi- 
tion); cf [llj and the references therein. At this level of (weak) regularity, 
the spacetime may contain impulsive gravitational waves (propagating at 
the speed of light) as well as shock waves (propagating at about the speed 
of sound in the matter). While shock waves occur only when the Einstein 
equations are coupled to the Euler equations of compressible fluids, the weak 
regularity class is relevant even in vacuum spacetimes and allow for singular 
nonlinear wave patterns not observed in the regular class and yet relevant 
from the physical standpoint. 

In recent work by the author and collaborators pQ H4L I11L I15L I12L 113] . 
several classes of spacetimes were constructed by solving the initial value 
problem for the Einstein equations when the initial data set has weak regu- 
larity and enjoys certain symmetry assumptions. This work encompasses es- 
pecially spacetimes admitting two Killing fields and satisfying the Einstein- 
Euler equations. 

All classical work on spacetimes with two Killing fields requires that 
data (and solutions) are sufficiently regular. When weak regularity, only, is 
assumed, many classical arguments and estimates are no longer valid, and 
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new arguments of proof must be developed that are robust enough to estab- 
lish the existence of weakly regular spacetimes and determine their global 
geometric structure. In comparison with classical results, by encompassing 
a larger class of initial data sets and, therefore, allowing for propagating 
waves of possibly singular nature, one need to go significantly beyond the 
techniques of analysis for regular spacetimes, and one provides additional 
rigorous validation of the underlying physical theory (including Penrose con- 
jecture) within a much broader context. 

An outline of this review is as follows. In Section [21 we consider general 
manifolds that need not admit symmetries, and we present a general def- 
inition of curvature for weakly regular metrics. In Section [3j we focus on 
spacetimes with Gowdy symmetry when the matter is an isothermal perfect 
fluid, and we present a global foliation of the future development of an initial 
data set with weak regularity. Finally, in Section [H we consider the special 
class of plane-symmetric spacetimes when the matter is a null fluid, and 
we establish a version of Penrose's conjecture for the broad class of weakly 
regular spacetimes. 

2. Definition and properties of manifolds with weak regularity 

Distributions on manifolds. We begin by presenting the definition of 
the curvature introduced in [11] . which is at the basis of our investigation. 
Our first task is to define distributions on an m-dimensional differential 
manifold M (with smooth topological structure). The definition is motivated 
by the canonical embedding / G L\ oc (M) i-> / G V(M), determined by 

<f,0J^v>,V-= [ foJ, w G DA m (M), 
Jm 

where T>A m (M) denotes the space of all (compactly supported) m-form fields 
on M (or densities, or by abuse of notation, volume forms). Precisely, the 
space of scalar distributions T>'{M) is defined as the continuous dual of 
T>A m {M). Analogously, the space of distribution densities V ' A m {M) is the 
continuous dual of the space T>(M) of all compactly supported functions on 
M. 

Given a smooth vector field X and an m-form field uj, we denote by 
ix^ the associated interior product. Given any open subset with regular 
boundary ft C M and using the identity Cx = dix + ixd together with 
Stokes formula, we obtain 

f (X/)oj = - [ f£ x u+ [ fi x u. 
Jn Jn Jan 

This observation motivates the following definition of the action Xf of a 

smooth vector field on a scalar distribution: 

(2) -< Xf, uj y v , !V := - -< f, C x u w G VA m (M). 

The notion of tensor distributions is defined as follows. Denote by 
£g(M) := C°°7f (M) the space of all smooth (p,g)-tensor fields on M, and 
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define the space VTq(M) of all (p, q)-tensor distributions as the space of all 
C°° (M)-multi- linear maps 

A : ZtiM) x ... x 1q(M) x T?(M) x . . . x T?(M) -> P'(M). 

V v ' V v ' 

q times p times 

For instance, the embedding L\ oc Tq(M) C V'Tq(M) is canonically defined 
by 

, LO >-T>',T> 

:= / A(X {1) ,...,X (q) ,e^,...,e^)co 

JM 

for all uj G 2?A m (M), , . . . , G Tj(M) and 0W, . . . ,0(p) € T?(M). 

Connections on manifolds. Adistributional connection V : Xg(M) x 
Tq(M) — >■ P'Tq (M), by definition, satisfies the standard linearity and Leib- 
nitz properties for all smooth fields. However, this class is often too broad 
in the applications, and more regularity may be assumed. 

An Lf oc connection, by definition, has V xY G Lf oc T^(M) for all X, Y G 
To(M), and its canonical extension: V : Xj(M) x Lf oc T^(M) ->• P'T^M) 
is defined by 

(v x y #} = x«y, 0)) - <y, V X 0> in V(M), 

XgTJ(M), YeLlX(M), 9gVT?(M). 

In this class of connections, we define the distributional Riemann tensor 

Rm : ll (M) x Tj(M) x Tj(M) -> P / T 1 (M), 

as an element of the dual space Hr} by 

(Rm(I, Y)Z, 9) =X(V Y Z, 9) - Y(V X Z, 9) 

- {V Y Z, V x 9) + (V X Z, V Y 9) - (V [XX ]Z, 9) 

for all 9 G $?(M), X,y,Z G Sq(M). By taking the trace in the usual 
manner, we also define the distributional Ricci tensor Ric in Hr}. 

Sequences of connections. One can establish the following stability 
property: if V^ n ^ is a sequence of Lf oc connections on M, converging in Lj 2 oc 
to some connection V^ 00 -*, i.e. 

vj?Y -+ V^V in Ll c 

for all X, Y G Tq(M), then the distributional Riemann and Ricci tensors 
Rm'"' and Ric^ of the connections converge to the distributional curva- 
ture tensors of the limit, i.e. 

RraW ->. Rm' 00 ' Ric^ -»• Ric^ 

in the distribution sense. 
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Distributional Levi-Cevita connections. Consider next the general 
class of distributional metrics. Observe that the compatibility condition 
Vg = between a metric and its associated Levi-Cevita connection must be 
handled with care, since non-smooth connections do not act on non-smooth 
tensors. 

The distributional Levi-Cevita connection of a distributional metric g is 
the operator V b : (X,Y) G %(M) x Tj(M) ^ V X Y G V'T?(M), defined 
by the "dual" Koszul formula 

(V X Y, Z) := \ (x(g(Y, Z)) + Y(g(X, Z)) - Z(g(X, Y)) 

(3) 2 V , 

- g(X, [Y, Z]) - g(Y, [X, Z}) + g(Z, [X, Y})) . 

It follows that in the weak sense 

V X Y -V Y X -[X,Yf = 0, 

X(g(Y, Z)) - (V X Y, Z) - (Y, V X Z) = 

for all X,Y,Z G %{M). 

One can check the following stability property under distributional con- 
vergence. Let be a sequence of distributional metrics converging in V 

to some limiting metric g^°°^. Then, the distributional connection as- 
sociated with g^ converges in V to the connection associated with 
g(°°\ i.e., for all X,Y,Z G Xj(M) 

<V b 5V, Z) -> <V b 5°V, Z> in P'(M). 

Distributional curvature of metrics. In practice, further regularity 
is often assumed and, especially, the class of H l regular metrics plays an 
important role. 

If V is the Levi-Cevita connection of a metric g of class H^ oc n L^ c which 
is also uniformly non-degenerate (for some c > 0) 

(4) |det(<7)|>c, 

then it is of class Lf oc and its Riemann and Ricci curvature tensors are well- 
defined as distributions. Moreover, its scalar curvature R is also well-defined 
as a distribution: 

R = g^mc(E {a) ,E m ), 
where Ei a \ is a local frame. 

Namely, since H^ oc r\L^ c is an algebra and g is uniformly non-degenerate, 
it follows that g al3 G H^ oc n L~ C (M). On the other hand, one can write 

Ric(£ (a) ,£ (/3) ) = (E^,Rm(E (a) ,E (a) )E W) ) 

= E (a) , V £(ct) > ) - E {a) , V E(a) E W > ) 

" <V B(a) £<"\ V S(CT) %)> + (VE w Et"\V Ela) E w ) 
- (E^,V [E(a)jE(a)] E w ). 
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Since the last three terms belong to L\ oc (M), we only need to define the 
product of g a P and the distributions 

E [a) «£? W , Vu w E W ) ) , ((E^ , V B(a) %) > ) . 

This is done by letting (for the first term, say, since the second term is 
similar) 

9 ^E {a) ((E^,V E{a) E m )) 

:= E {a) (g^(E^,V E(a) E (p) )) - (^(^.V^), 
which, clearly, is well-defined as a distribution. 

Sequences of metrics. One can check the following property of stabil- 
ity under strong sconvergence. Let g^ be a sequence of H} metric tensors 
converging locally in H} to some limiting metric g^ 00 " 1 , and assume that the 
inverse metrics gT-, converge locally in Lf£ c to gT^y- 

• The Levi-Cevita connections associated with g( n J are of class 
Lf oc and converge in Lf oc to the connection V*- 00 - 1 of the limit g^°°\ 
that is, for all X, Y € 1^(M), 

Y -»■ V^° ) y strongly in L? oc . 

• The distributional Riemann, Ricci, and scalar curvature tensors 
RmW, Ric (n) , of the connections V (n) converge in T>' to the 
limiting curvature tensors Rm' m ', Ric* 00 * 1 , R^ of V^°°^ respec- 
tively. 

Hence, the curvature of a spacetime may be defined in the sense of dis- 
tributions when the metric has weak regularity, only. We have advocated 
here the use of fully geometric definitions, while coordinate dependent def- 
initions are usually adopted in the literature; see Geroch and Traschen [9], 
Mars and Senovilla |17j . and Lichnerowicz |16j . Other approaches are de- 
veloped for certain classes of spacetimes by Kunzinger, Steinbauer, Vickers, 
Balasin, Aichelburg, Heinzle, etc. 

A framework to deal with manifolds with low regularity was also intro- 
duced by Lott, Villani, Sturm, and allows to define Ricci lower bounds via 
the theory of optimal transport. Both lower and upper curvature bounds 
can also be defined in the theory of Alexandrov spaces. 

Discontinuous connections. Consider connections having a jump dis- 
continuity across a given smooth hypersurface H C M. A connection 
that is discontinuous across H is, by definition, a connection V of class 
tfJM) n Wl£{M±) having a jump discontinuity across H, with 

M = M~UM + , M~C\M + = U. 
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Clearly, according to our definitions above, the corresponding operators V ± 
determined from V by restriction to each M ± have well-defined, distribu- 
tional Riemann and Ricci curvatures Rm ± and Ric^. In fact, these cur- 
vature tensors belong to L 1 1 oc (M ± ), at least, and actually to V\ oc {M^) if 
p > m/2. 

Our purpose now is to compute the distributional curvature of the con- 
nection V, which we expect to contain Dirac mass singularities along H. 

First of all, given a smooth hypersurface H, the Dirac measure supported 
by U C M is defined as the 1-form distribution 8 H & V"lf(M) such that 

X e%l(M) ^ (d H ,X) eV(M), 
<{8 H ,X),u^v>,v= [ ixoo, oj£VA m (M). 

JH 

We write \A\u '■= A + — A~ for the jump of a tensor field A across W, and 
we write the "regular part" as 



A reg := 



A+ in M+, 
A~ in M~. 



For instance, the distributional derivative W of a vector field (smooth 
in and discontinuous across T~L) is 

(5) VxV := (VxVY^ + [V] H (S n ,X) 

for all X G 1q(M). The term (8 H ,X) depends on X only via its restriction 
to H, so 8% can be applied to vector fields only defined on %. If X is a 
vector field tangent to the hypersurface (X G Tq(%)), then = 0. 

Jump relations. We can now derive jump relations associated with 
a discontinuous connection. Choose an adapted frame so that E^, i = 
1, . . . , m — 1 is a local frame on the hypersurface H, while , a = 1, . . . , m: 
is a local frame on M, and let E^: be the associated dual frame of 1-form 
fields. Then, the distributional Riemann curvature takes the form 

Rm(X, Y)Z = {Rm(X,Y)Z)™s + [V Y Z] n (8 n ,X) 

-[VxZ]h(Sh,Y) 

for all X, Y, Z € 1q(M), while the distributional Ricci curvature reads 

Ric(X,y) = (Ric(X,Y)Y^ + j V E(a) Y)] H (8 n ,X) 

-[(E( m \v x Y)] n (8n,E (m) ) 

for all X,Y e %1{M). 

From these formulas, we can deduce some properties of the jumps. Con- 
sider a connection V that is discontinuous across a hypersurface Ti. 

• The singular part of the Riemann tensor vanishes if and only if the 
connection V is continuous across %. 
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• The singular part of the Ricci tensor vanishes if and only if the 
components 

^E^ m \V E j)X^ and (E^\ V E (j) X) are continuous across % 
for all vector field X G lJ(M). 

• Suppose V is the L 2 Levi-Cevita connection of a uniformly non- 
degenerate metric g and suppose that g is on each side of 
some hypersurface H, then 

R := R re9 + [(g^E® - g^E^\ V E{j) E w )] n (dn, E {m) ). 
where E^' is a frame adapted to H. 

3. A global foliation for Einstein-Euler spacetimes with 
Gowdy-symmetry on T 3 

Spacetimes with Gowdy symmetry. This section is based on the 
paper |12j in collaboration with A.D. Rendall. We impose that the space- 
time (M,g) under consideration has Gowdy symmetry, which is a classical 
assumption made to study inhomogeneous cosmology with a "big bang" 
or "big crunch". As far as the matter model is concerned, we consider a 
compressible, isothermal perfect fluid described by the energy-momentum 
tensor 

(6) T a/3 = {fi + p)u a u 13 +pg a(3 

with future oriented, timelike velocity vector u a , normalized so that 

g a pu a vP = -1, 

and with mass-energy density /i > and pressure p = c 2 s \x. Here, c s € (0, 1) 
does not exceed the light speed normalized to be 1. Here, all greek indices 
a, /3, . . . describe 0, 1, 2, 3. 

Our main objective here is to construct a class of globally foliated space- 
times, obtained by solving the initial value problem when initial data have 
low regularity. Under our construction, gravitational waves (possibly impul- 
sive ones) at the speed of light can propagate in such spacetimes, and shock 
waves can also propagate at (about) the speed of sound. 

Gowdy symmetry on T 3 imposes that the (3 + l)-dimensional Lorentzian 
manifold (M,g) admits the Lie group T 2 as an isometry group acting on 
the torus T 3 : it is generated by two (linearly independent, orthogonally 
transitive) spacelike vector fields X, Y satisfying, by definition, the Killing 
condition 

C x g = Cyg = 0, 

the commutation condition 

[X,Y]=0, 

the spacelike property 

g(X,X),g(Y,Y)>0. 
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In addition, by definition, we impose the property of orthogonal transitivity, 
in the sense that the distribution of 2-planes of covectors 

Vect(X,F) ± 

is Frobenius integrable. 

In the vacuum, the above conditions determine exactly the so-called 
Gowdy spacetimes. Note that the above property of integrability can be 
analyzed as follows. By introducing two vectors Z, T orthogonal to X, Y, we 
can ensure that (using a reduction to dimension 2 and after normalization) 

e(X,Y,-,-f = Z®T -T®Z, 
where e is the canonical volume form, or equivalently 

e aPrf5 X^Y s = Z a T? -T a Z p . 

Then, the distribution of covectors g(X, •), g(Y, •) is Frobenius integrable if 
and only if Z, T, [Z,T] are linearly dependent, that is, e(Z,T, [Z,T], •) = 
or 

e aM5 Z a TP[Z,TY< = 0. 

After some calculations, this condition is found to be equivalent to vanishing 
"twist constants" (as identified by Chrusciel [6]) 

(e(X, Y, •), VX) g = (e(X, Y, ; -),VY) g = 0, 

i.e. 

e a ^ s X a Y^VX 5 = e a ^ s X a Y^'Y 5 = 0. 

Note that the vacuum Einstein equations imply that the above (twist) quan- 
tities are constants, while this is an assumption in matter spacetimes. 

Main existence result. A Gowdy symmetric initial data set on T 3 

(H,g,k, p,j) 

is determined by a Riemannian 3-manifold (H,g) together with a symmetric 
tensor field k (playing the role of the second fundamental form of the space- 
time to be constructed), in which ~p denotes the energy density of the fluid 
(a scalar field) and j its momentum (a vector field, essentially equivalent 
to the velocity field) defined on H. In addition, the so-called Einstein con- 
straint equations must hold, and these data must be invariant under Gowdy 
symmetry. These initial data is assumed to be solely weakly regular. 

The class of weakly regular spacetimes with Gowdy symmetry on T 3 is 
defined by the following conditions: 

• The metric belong to the Sobolev space H (£) on every spacelike 
slices £ C M. 

• The mass-energy density p > (scalar field) and the momentum 
j (vector field) belongs to the Lebesgue space ^ 1 (S) of integrable 
functions, for every spacelike slice X. (Prescribing p,j is equivalent 
to prescribing p,u.) 
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• The Einstein equations ([I]) are satisfied in the distribution sense 
and, in particular, the curvature is well-defined as a distribution in 
-ff -1 (£) for each slice S. 

• Furthermore, the entropy inequalities associated with the Euler 
equations are satisfied. 

Theorem 3.1 (Einstein-Euler spacetimes with Gowdy symmetry on T 3 ). 
Let (g, k,p, j) be a weakly regular, Gowdy symmetric, initial data set on T 3 
for the Einstein-Euler equations, and assume that these initial data have con- 
stant area, and are everywhere expanding or everywhere contracting. Then, 
there exists a weakly regular, Gowdy symmetric spacetime (M,g,p,j) satis- 
fying the Einstein-Euler equations in the distributional sense, so that the fol- 
lowing properties hold. The manifold (M, g, p, j) is (up to diffeomorphisms) 
a Gowdy-symmetric future development of (g,k,p,j) which is globally cov- 
ered by a single chart of coordinates t and (6,x,y) € T 3 , with 

{(i, 6) / < Co < t < 00} x T 2 , expanding case, 

M = I 

{(t,0) j Co < t < c\ < 0} x T 2 , contracting case. 

Here, c\ € (co,0] is a constant, and the time variable is chosen to coincide 
with the area of the surface of symmetry in the expanding case, and with 
minus this area in the contracting case. 

Further comments and perspectives. Many earlier works on reg- 
ular spacetimes with Gowdy symmetry are available, both, in the vacuum 
case (Berger, Chrusciel, Isenberg, Moncrief, Weaver) and when the Einstein 
equations are coupled to a kinetic matter models (Rendall, Andreasson, 
Dafermos); cf. [18J and [12J for references. For an analysis of the global 
geometry of T 2 -symmetric vacuum spacetimes with weak regularity, we re- 
fer to LeFloch and Smulevici [13J, and for an earlier result on compressible 
fluid models (a "small time" result for spacetimes with bounded variation) 
to LeFloch and Stewart [15J (see also [lj). 

Note that the time- function constructed above is geometric in nature and 
is based on the Gowdy symmetry assumption (areal coordinates). It would 
be interesting to investigate whether a foliation by constant mean curvature 
slices also exists. Even more important is to investigate the structure of 
the boundary of the future development of the constructed spacetimes. In 
general, the constant c\ need not vanish, since explicit counter-examples 
for the Einstein-Euler equations are available, although such examples are 
probably unstable. 

On the other hand, recall that spacetimes with bounded variation were 
constructed by Christodoulou in his work [3], |4j that settled positively the 
weak version of Penrose's cosmic censorship conjecture in the context of 
spherically symmetric spacetimes and for scalar fields. Recall also that 
Groah and Temple |10j established a local-in-time existence result for spher- 
ically symmetric matter spacetimes. In such spacetimes, no gravitational 
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waves are permitted and the matter equations are coupled with a differen- 
tial equation accounting for (non-evolutive) geometrical features. 

Following Gowdy and Geroch, one considers the quotient manifold M := 
M/T 2 . Since the two Killing fields (X,Y) =: (Xi, X2) are linearly indepen- 
dent and spacelike, the area of the orbit of orbits of the T 2 symmetry group 
satisfies 

(7) R 2 := det (A a6 ) > 0, X ab := g(X a , X b ). 

Due to the Killing field property, R > is a constant on each orbit. 

In addition, the associated Lorentzian metric defined on M can be ex- 
pressed as h := g — X ab X a <£> X b . and one can also introduce the projection 
operator: 

h v „ u \ ab Y Y u 

n fi ■— 9fi — A -A-afi^b ) 

which serves to project the Einstein equations on the quotient manifold. 

This leads us to introduce a set of local coordinates (t, 9, x, y) such that 
x, y span the orbits of symmetry, and the metric coefficients depend on (t, 9), 
only, and are periodic in 8. A key property that allows to choose the area 
as a time-function is that the gradient Vi? is timelike — a fact that remains 
true for the Euler equations and even in the weak regularity class under 
consideration. 

More precisely, in suitably conformal coordinates (r,9,x,y) one can ar- 
range that the metric reads 

(8) g = e 2 ^-U) (-dr 2 + d9 2 ) + e 2U {dx + A dy) 2 + e~ 2U R 2 dy 2 , 

where 77, U, A, R with R > depend on r, 9, and the area R is one of the un- 
knowns. The Einstein equations reduce to a coupled system of eight partial 
differential equations, consisting of four evolution equations (second-order, 
semi-linear wave equations for U, A, r], R) and two constraint equations (non- 
linear differential equations for Vi?), which also imply the Euler equations 
(nonlinear hyperbolic equations for p,j, or fi,u). 
In area] coordinates (t,9,x,y) one has 

(9) g = e 2 ^ {-dt 2 + a~ 1 d9 2 ) + e 2U (dx + A dy) 2 + e~ 2U t 2 dy 2 , 

in which U, A, is, a depend on t, 9, and the time t = R coincides with the 
area of the orbits of symmetry. The Einstein equations reduce to a cou- 
pled system of six partial differential equations, specifically three evolution 
equations (second-order, nonlinear wave equations for U, A, v) and three con- 
straint equations (including nonlinear differential equations for a), while the 
Euler equations again form a system of two nonlinear hyperbolic equations. 

4. The characteristic initial value problem for plane symmetric 
spacetimes with weak regularity 

Objective in this section. Based on the paper |15j in collaboration 
with J.M. Stewart, we consider the class of spacetimes that satisfy the "po- 
larized" Gowdy symmetry property (hypersurface orthogonal property), and 
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we investigate the Einstein equations when the matter is a null perfect fluid 
(cf. (HO]), below). 

We are interested here in the characteristic initial value problem when 
data are prescribed on two null hyper surf aces, denoted by N UQ and M Vo , 
intersecting along a 2-plane V uo , vo - Like in the previous section, spacetimes 
with weak regularity are constructed by solving the initial value problem. 
Our set-up and assumptions allow us to tackle a new issue, beyond the 
sole existence of the spacetimes, that is, to establish a suitable version of 
Penrose's strong censorship conjecture in this class. 

Recall that the characteristic initial value problem addressed here for 
matter spacetimes with symmetry and weak regularity, was treated earlier 
within the class of sufficiently regular vacuum spacetimes by Friedrich [8] . 
Stewart and Friedrich [20J , Rendall |19j , Christodoulou [5] , Dossa and Tad- 
mon [7], and Choquet-Bruhat, Chrusciel, and Martm-Garcfa [2]. 

Matter model of interest. The matter model under consideration is 
a null, irrotational fluid, whose sound speed, by definition, coincides with 
the speed of light (normalized to unit). In this context, propagating dis- 
continuities in the fluid variables may occur only along null hypersurfaces. 
More precisely, a perfect fluid with pressure p = [i equal to its mass-energy 
density [i (null fluid) has energy-momentum tensor 

(10) T a ? = 2fiu a u p + fig alB . 

As usual, u is the velocity vector, normalized so that u a u a = — 1. 

The (second) contracted Bianchi identity allow us to derive the Euler 
equations 

V Q T a/3 = 0, 

which take the form 

(11) (u a V a fi) v? + n (V a u a ) v? + fiu a V a vP - -V^// = 0. 
Multiplying (llip by up we obtain the scalar equation 

2V a w u a - 2w V a u a + 2w u a upV a u p - V ' a w u a = 0, 

which, in view of upW a u^ = 0, simplifies 

u a V a w + 2wV a u a = 0, 

or equivalently, assuming that the density is bounded away from zero and 
setting S = i logiu, we obtain 

n a V Q S + V a u a = 0. 

Next, multiply by the projection operator iJg 7 = — ugii 7 and obtain 
the vector equation 

H a ~<V a w - 2w u a V a u< = 0, 

or equivalently 

iT 7 £ a - u a V a u^ = 0. 
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In addition, we impose that the flow is irrotational, that is, we assume 
the existence of a (scalar) potential ip with timelike gradient 

v«v vfy < 0, u a = 1 . 

After normalization (by replacing tp with F(ip)), the Euler equations imply 
the so-called Bernouilli's equation (in characteristic variables, cf. below) 

which yields the energy density in terms of the potential. We arrive at a 
single matter equation for the scalar field ip. In characteristic coordinates, 
this is a wave equation for the velocity potential 

ipuv + b v ipu + Ki>v = o. 

Reduced form of the field equations. Each Killing field is assumed 
to be hypersurface orthogonal ("plane symmetry") and we arrange the co- 
ordinates so that the metric reads 

g = e ^ (-dt 2 + dx 2 ) + e 2b (e 2c dy 2 + e~ 2c dz 2 ) 
= -e 2a dudv + e 2b (e 2c dy 2 + e~ 2c dz 2 ), 

where the coefficients a, b, c depend upon the characteristic variables 

u = t — X, V = t + X, 

only. The relevant components of the Einstein tensor are 

G 00 = 2 (-2 a u b u + b uu + b 2 u + c 2 u ) , 

Gq\ = 2 (—b uv — 2 b u b v ) , 

Gn = 2 (-2 a v b v + b vv + b 2 v + c 2 v ) , 

G22 = 4e^ a+/0+/c ( + 
G33 = Ae- 2a+2b - 2c (a uv + b uv + b u b v o u c v + b 

We use here one of our gauge freedom, and solve a decoupled equation 
for b, allowing us to choose coordinates so that 

e 2b = \\u + v\, 

and the region of interest is 

{u + v < 0}, 

the hypersurface u + v = being a (possibly physical or coordinate) singu- 
larity. 
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The corresponding evolution equations are two singular wave equations 
of Euler-Poisson-Darboux type 

(12) 2{ \ +v) 

C-uv ~\~ 777 ; r ( C-u ~\~ C v ) = , 
2{u + v) 

while the constraint equations read 

„2 i 1 „,.2 



(13) 



a, = (cg + ^)(n + ^)- 4( ^ 1 + ^ ) . 



Characteristic initial value problem. 

We can establish a global and fully geometric result, based on prescribing 
initial data on two null hypersurfaces intersecting along a two-plane. We 
seek for the global structure of the future development of a given initial data 
set. We consider arbitrary null coordinates, denoted below by (U, V) and, 
throughout, we restrict attention to plane symmetric data and spacetimes. 

An initial data with weak regularity consists of the following prescribed 
data. Let 

(M, e^dUdydz) , (Jf, e^dVdydz) 
be two plane symmetric 3-manifolds (endowed with volume forms) with 
boundaries identified along a two-plane V and parametrized for some (Uo, Vo) 
as 

M:={U>U }, Jf:={V>V }, V := {U = U , V = V } . 

(1) Suppose that a, a belong to the Sobolev space W 1 ' 1 , i.e. the inte- 
grals 



/ 

J N 



(\a\ + \dua\)e^dU [ (\a\ + \d v a\) e a dV 



!N_ JJT 
are finite, and are normalized so that a\j> = a\p = 0. 
(2) Let $oo an d $4> $22 be (plane-symmetric) functions defined on 
the hypersurfaces M_ and AT, respectively, with < $ 00 G L 1 (A/') 
and < $22 € L l (N), i.e. the integrals 

/ $oo e-dU, [ $ 22 e"dF 

are finite, and that _3> G W' 1 ' 2 ^) and $ 4 G W -1 ' 2 ^), i.e. * = 
du^o ] and "^4 = dv^4 ] with 



/ 



\*$\ 2 e£dU, [ <£%\ 2 e a dV. 



/AT JaT 
(3) Finally, one also prescribes the connection NP scalars p$, uq, \q, fiQ 
on V. 
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We then have the following existence and qualitative property result, for 
which we refer the reader to [15] about the terminology used here. 

Theorem 4.1 (Global causal structure of plane symmetric matter space- 
times). Consider an initial data set with weak regularity determined by 
(M, e 2 -) and (If, e°) , a plane (V, po, do, Ao, po), and prescribed Ricci and 
Weyl NP scalars (* !^oo) and (*4>$22). 

(1) Then, there exists a W 1,2 regular spacetime (A4,g) determined by 
metric coefficients a, b, c and matter potential ip which is a future develop- 
ment of the initial data set satisfying the Einstein equations ([1]) for self- 
gravitating, irrotational fluids with the initial conditions 

(p,o-,X,p) = (po,cr ,X Q ,Ho) onV, 
[a, ^q, $00) = («, 3£o, ^00) 071 the null hypersurface Af , 

and 

(a, ^4, $22) = (s, x I / 4, ^22) on the null hypersurface J7. 

(2) The constructed development of the initial data has past boundary 

{Ho>U > U ;V = Vo) U {U = U ;V > V > V } c AAuAT 

and, for generic initial data (in the sense defined in (3), below) the cur- 
vature blows up to (and makes no sense, even as a distribution) as one 
approaches its (W 1 ' 2 regular) future boundary 

Bo := {F(U) + G(V) = 0} 

for some functions F, G, so that the spacetime is inextendible beyond Bo 
within the class of W 1 ' 2 regular spacetimes . 

(3) The above result holds for generic initial data, in the sense that 
arbitrary data can always be perturbed in the natural (energy-type) norm 
so that the perturbed initial data do generate a singular spacetime whose 
curvature blows-up on Bo- 

The coefficients e— and e a , modulo a conformal transformation, could 
be chosen to be identically 1 on the initial hypersurface, so that two main 
degrees of freedom remain on each of the two initial hypersurfaces. Theo- 
rem l4,ll can be seen as a statement of Penrose strong's censorship conjecture 
for plane symmetric spacetimes with weak regularity. 
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